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In the last quarter of a century, a unified characterization of the spectral evolution of low-mass 
X-ray binaries, both containing a neutron star and a black hole, was possible. In this context, the 
notion of source states characterizing the X-ray emission from black-hole binaries and neutron- 
II star low-mass X-ray binaries revealed to be a very useful tool to disentangle the complex spectral 

O and aperiodic phenomenology displayed by those classes of accreting objects. Be/X-ray binaries 

constitute another major class of transient accreting binaries, for which very little work has been 
Ci done on the correlated timing and spectral variability. Especially, no definition of source states 

exists for this class, in spite of their highly variable X-ray emission. When active, Be/X-ray bina- 
I ries are among the brightest objects in the X-ray sky and are characterized by dramatic variability 

K*' in brightness on timescales ranging from seconds to years. It is then worth it to ask whether a 

definition of spectral states is possible for these systems. 
\G In this work, we try to address such a question, investigating whether accreting X-ray pulsars 

display source states and characterizing those states through their spectral properties. Our results 
show that Be/X-ray pulsars trace two different branches in their hardness-intensity diagram; the 
horizontal branch, a low-intensity state, and the diagonal branch, a high-intensity state that only 
I appears when the X-ray luminosity exceeds a critical limit. We propose that the two branches are 

. 5_| the phenomenological signature of two different accretion modes - in agreement with recently 

proposed models - depending on whether the luminosity of the source is above or below a critical 
value. 
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1. Introduction 

In the last twenty-five years, a unified description of tlie spectral evolution of low-mass X-ray 
binaries (LMXB) and black-hole binaries (BHB) was made possible [1-4]. In such systematic 
studies, color-color diagrams (CDs) and hardness-intensity diagrams (HIDs) emerged as powerful 
tools to characterize those systems, since it was observed that, both for persistent and transient 
systems, the source emission evolves smoothly along the diagrams, describing a pattern. Similarly 
to color- magnitude diagrams in the optical (although with much faster time scales), CDs and HIDs 
permitted to unveil and address questions that could not be managed otherwise, by introducing the 
crucial notion of spectral states. A state is defined by the specific and well-defined position that 
the sources occupy in their CDs/HIDs, and is associated to the appearance (or disappearance) of a 
characteristic spectral or timing component in the corresponding energy spectra and power spectral 
density. 

Within high-mass X-ray binaries (HMXBs), Be/X-ray binaries (Be/XRBs) represent the major 
subclass, and one of the most extreme case of X-ray variability on several time scales. These 
systems consist in a strongly magnetized neutron star orbiting a massive Be star. The observed 
X-ray emission is powered by large amount of matter accreted onto the compact object, during two 
different classes of events: 

- type-I (or normal) outbursts, lasting a few days, with luminosities ~10^^- 10^^ erg s^^ occur 
in series separated by the orbital period, generally close to the time of periastron passage of 
the neutron star; 

- type-II (or giant) outbursts, rarer and unpredictable, last several weeks, showing 2-3 orders 
of magnitude variations in X-ray intensity (Lx > 10^^ erg s^^). 

This work focuses on the second type of X-ray outbursts, with the purpose of investigating 
the spectral variability in Be/XRBs as a class, and the final objective of defining and characterizing 
spectral states. A first attempt to define spectral states through color analysis was carried out by 
[5, 6]. The present work builds up on those results, extends them to include all the Be/XRBs de- 
tected by RXTE during a type-II outburst, presents a finer sampling resolution in spectral analysis, 
and provides an interpretation of the identified states within the recently proposed scenario that 
describes accretion in Be/XRBs. 

2. Analysis 

We analyzed archived data from the Rossi X-ray Timing Explorer (RXTE). The hardness, or 
soft color (SC), was obtained directly from the Proportional Counter Array (PCA) count rates in 
the 7-10 keV / 4-7 keV energy bands. Energy spectra were extracted, one spectrum per pointing, 
for the PCA and the High-energy Timing Experiment (HEXTE), and simultaneously fitted. The 
lack of adequate theoretical continuum models for accreting neutron stars implies the use of em- 
pirical models to fit the observations. To compare the results from all sources consistently, we 
used the same spectral components for all the system analyzed. The best-fit model consisted in a 
combination of photoelectric absorption and a power-law with high-energy exponential cutoff to 
fit the continuum. The discrete components were accounted for with an emission Gaussian line 
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Source 


Time range 


Num. 


Exp. 


Max. Lx 


Ecyc 


7 


name 


(MJD) 


obs 


time (ks) 


(ergs ') 


(keV) 




4U Ul ID+OJ 


J jzj4. 1— ji JU4.0 


36 


96.4 


1 /I V/ 1 n38 
1 .4X 10 


12 


1 .5 ± 0.4 


KS 1947+300 


51874.2-52078.0 


86 


144.9 


7.1x10" 




1.1±0.4 


EXO 2030+375 


53914.9-54069.9 


154 


140.1 


1.6x10^'* 


11 


l.liO.l 


V 0332+53 


53340.3-53376.6 


80 


174.5 


3.4x103** 


25 


1.6±0.5 


lA 1118-616 


54838.3-54865.0 


26 


10.2 


2.9x10^' 


60 


1.2±0.2 


XTEJ0658-073 


52932.7-53056.4 


76 


201.9 


3.6x10^' 


33 


1.1±0.3 


Swift J1656.6-5156 


53723.9-54269.9 


260 


313.0 


5.2x10" 




0.9±0.2 


lA 0535+262 


55169.4-55210.5 


27 


126.7 


6.8x10" 


46 


0.7 ±0.2 


GROJ1008-57 


54426.0-54457.2 


20 


27.5 


2.0x10" 




0.8±0.2 



Table 1: RXTE observations. 



profile, which fits the iron Fe Ka fluorescence, and Lorentzian absorption profiles, which fitted 
cyclotron resonance scattering features (CRSF). For three sources, 4U 0115+63, XTE J0658-073 
and 1 A 0535+262, additional residuals around 7-10 keV required the inclusion of a broad Gaussian 
emission line, a feature known as the "bump" model [7]. Table 1 shows the journal of observations, 
including the maximum X-ray luminosity reached at the outburst peak, and a measure of the central 
energy of the CRSF, when present, and the average x^ed ^'^^ spectral fitting. 

3. Evidence for source states in Be/XRBs 

In this section we present our results, which unveil the existence of spectral states in X-ray 
pulsars from both model-independent color analysis and spectral analysis. 

3.1 Hardness-intensity diagrams 

For all the sources analyzed in this work.. Fig. 1 shows the corresponding HID, a plot of the 
3-30 keV band count rate as a function of hardness, one point per RXTE observation. In the fig- 
ure, a color shade marks the four sources which display two branches in their HID, with a clear 
transition between them: 4U 0115+63, EXO 2030+375, V 0332+53, and KS 1947+300. These 
correspond to the brightest sources (see Table 1), although the luminosity does not seem to be the 
only parameter triggering the state transition. For consistency with [5, 6], we will keep the termi- 
nology horizontal branch (HB) and diagonal branch (DB) to designate the low-luminosity branch 
and the high-luminosity branch, respectively, although in the present work, as we are plotting the 
HID using a logarithmic representation, both branches appear as diagonal in the diagram. 

Due to the unpredictable nature of type-II outbursts, it is usually difficult to have a proper 
coverage of the very first phase of the rise. Nevertheless, for one system, KS 1947+300, data were 
collected during the full duration of the outburst. We shall thus take KS 1947+300 as a paradigm of 
the objects displaying the two branches, and use it as an example to describe the pattern described 
by the source during the outburst. 

Figure 2 shows the outburst profile and HID for KS 1947+300, where we marked with different 
colors/symbols the different phases of the outburst. When the outburst begins, as the flux increases, 
the SC increases, i.e. the emission from the source becomes harder and harder, and the system is 
moving rightward on the HID, depicting the HB (red circles). When the luminosity reaches a 
critical value, the source undergoes the transition, and makes a sudden turn in the HID, entering 
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Figure 1: Hardness-intensity diagrams for the nine systems studied in this work. 

the DB (black triangles). At this stage, as the flux increases, the SC decreases, i.e. X-ray emission 
becomes softer and softer, and the system moves leftward on the HID, until the peak luminosity is 
reached. Then, as flux starts decreasing, the hardness starts increasing again, and the system moves 
back on the DB, tracing the same pattern as during the rise (blue squares). 




51900 52000 52100 0.54 0.56 0.58 0.60 0.62 0.64 0.66 0.68 0.70 0.72 0.74 0.76 

MJD SC 



Figure 2: Light curve (left panel) and HID (right panel) for KS 1947+300 

The transition occurs again at about the same critical luminosity as in the rise, when the source 
leaves the DB to enter the HB, and the hardness decreases as the flux decreases (red squares), back 
to the starting point in the HID. The two branches are drawn in exactly the same way during the 
rise and the decay, depending only on the luminosity and not on the outburst rise/decay phase, i.e., 
no hysteresis is observed. 

3.2 Spectral analysis 

The main parameter representing the spectral continuum in the spectral continuum in X-ray 
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pulsars in the classical X-ray range is the photon index T. Figure 3 shows the relation between T 
and luminosity, where Ledd = 1.7 x 10^^ erg s^^ is the Eddington luminosity for a neutron star. 
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Figure 3: Photon index as a function of X-ray luminosity. 



The same four systems that show two branches in their HID, also display a bimodal relation 
of r with luminosity, passing from a negative correlation to a positive correlation, after a critical 
luminosity is reached, confirming the existence of two luminosity spectral states. Instead, only one 
state is shown by the remaining systems. In Fig. 3, we marked with filled/open circles the pointings 
corresponding to the HB/DB respectively, as defined by the HID. One can see that there is fully 
agreement between results from color and spectral analysis. Although correlations between the 
photon index and luminosity have been reported in the past, we report for the first time a sudden 
change in the correlation sign, depending on the position of the source in its HID. 



4. Spectral states as a signature of a bimodal accretion 

Spectral variability in X-ray pulsars is thought to arise close to the neutron star surface by 
means of bulk and thermal Comptonization in an accretion column. This would form in correspon- 
dence of one or both polar caps, where the accretion flow is channelled by the neutron star mag- 
netosphere [8]. A detailed description of the interaction between the highly magnetized plasma 
surrounding the pulsar and the accretion flow is still lacking, although recently, an attempt to phe- 
nomenologically model the vicinity of the neutron star surface was able to describe the bimodal 
con^elation observed between the CRSF energy and luminosity [9]. Those authors invoked two 
different regimes of accretion depending on a critical luminosity, Ldit. Below (but close to) Lcrit, 
the deceleration of the accreting flow at the neutron star surface would be mainly due to Coulomb 
interaction, while above that value, the main decelerating factor would be radiation pressure. In 
this work we found that X-ray pulsars display two spectral states, depending on the luminosity 
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they reach during an outburst. Only the brightest sources enter the DB. We measured the transi- 
tional luminosity in the range ~ 1^ x lO''^ erg s^\ depending on the source. We propose that the 
transitional luminosity correspond to the critical luminosity in [9], and that the two branches we 
observe are a signature of the two accretion regimes, being the HB the sub-critical state, and the 
DB the super-critical one. In [9], the authors propose that the critical luminosity is proportional to 
the neutron star magnetic field intensity, or equivalently, to the energy of the cyclotron line. This 
is consistent with what we found. In fact, KS 1947-1-300 and lA 0535H-262 reached almost the 
same luminosity (see Table 1), but only the first system enters the DB. This behavior would be 
conveniently explained by the dependence of Lent on the magnetic field, which for 1 A 0535-1-262 is 
much higher and would thus prevent the source to undergo the transition. 

5. Conclusions 

We have performed color and spectral analysis of nine Be/XRBs during a giant outburst. The 
evolution through the hardness-intensity diagram suggests that Be/XRBs undergo state transitions, 
displaying, if a critical luminosity is reached, two different source states. These results are con^ob- 
orated by spectral analysis. The two states we observe are a plausible evidence of the two accretion 
regimes proposed by [9]. For a more detailed analysis of the results presented here the reader is 
referred to [10] 
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